Direct numerical simulation of turbulent channel flow up to Reτ≈5200Re_\tau \approx 5200

A direct numerical simulation of incompressible channel flow at $Re_\tau$ = 5186 has been performed, and the flow exhibits a number of the characteristics of high Reynolds number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Karman constant $\kappa = 0.384 \pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in $k$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bi-modal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.Comment: Under consideration for publication in J. Fluid Mec

Characteristic eddy decomposition of turbulence in a channel

The proper orthogonal decomposition technique (Lumley's decomposition) is applied to the turbulent flow in a channel to extract coherent structures by decomposing the velocity field into characteristic eddies with random coefficients. In the homogeneous spatial directions, a generaliztion of the shot-noise expansion is used to determine the characteristic eddies. In this expansion, the Fourier coefficients of the characteristic eddy cannot be obtained from the second-order statistics. Three different techniques are used to determine the phases of these coefficients. They are based on: (1) the bispectrum, (2) a spatial compactness requirement, and (3) a functional continuity argument. Results from these three techniques are found to be similar in most respects. The implications of these techniques and the shot-noise expansion are discussed. The dominant eddy is found to contribute as much as 76 percent to the turbulent kinetic energy. In both 2D and 3D, the characteristic eddies consist of an ejection region straddled by streamwise vortices that leave the wall in the very short streamwise distance of about 100 wall units

Coherent structures in a simulated turbulent mixing layer

A direct numerical simulation of a plane turbulent mixing layer has been performed. The simulation was initialized using two turbulent velocity fields obtained from direct numerical simulation of a turbulent boundary layer at momentum thickness Reynolds number 300 (Spalart, 1988). The mixing layer is allowed to evolve long enough for self-similar linear growth to occur, with the visual thickness Reynolds number reaching 14,000. The simulated flow is examined for evidence of the coherent structures expected in a mixing layer (rollers and rib vortices). Before the onset of self-similar growth, such structures are present with properties similar to the corresponding laminar or transitional structures. In the self-similar growth regime, however, only the rollers are present with no indication of rib vortices and no indication of conventional pairing. This results in a reduction of mixing and layer growth

The three-dimensional evolution of a plane mixing layer. Part 1: The Kelvin-Helmholtz roll-up

The Kelvin Helmholtz roll up of three dimensional, temporally evolving, plane mixing layers were simulated numerically. All simulations were begun from a few low wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity profile. The spanwise disturbance wavelength was taken to be less than or equal to the streamwise wavelength associated with the Kelvin Helmholtz roll up. A standard set of clean structures develop in most of the simulations. The spanwise vorticity rolls up into a corrugated spanwise roller, with vortex stretching creating strong spanwise vorticity in a cup shaped region at the vends of the roller. Predominantly streamwise rib vortices develop in the braid region between the rollers. For sufficiently strong initial three dimensional disturbances, these ribs collapse into compact axisymmetric vortices. The rib vortex lines connect to neighboring ribs and are kinked in the opposite direction of the roller vortex lines. Because of this, these two sets of vortex lines remain distinct. For certain initial conditions, persistent ribs do not develop. In such cases the development of significant three dimensionality is delayed. When the initial three dimensional disturbance energy is about equal to, or less than, the two dimensional fundamental disturbance energy, the evolution of the three dimensional disturbance is nearly linear (with respect to the mean and the two dimensional disturbances), at least until the first Kelvin Helmholtz roll up is completed

The three-dimensional evolution of a plane mixing layer. Part 2: Pairing and transition to turbulence

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings was simulated numerically. Initial conditions for all simulations consisted of a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found both to inhibit the growth of infinitesimal three-dimensional disturbances and to trigger the transition to turbulence in highly three dimensional flows. The mechanisms responsible for the growth of three-dimensionality as well as the initial phases of the transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of span wise vorticity, which undergo a secondary roll up. Transition also produces an increase in the degree of scalar mixing, in agreement with experimental observations of mixing transition. Simulations were also conducted to investigate changes in span wise length scale that may occur in response to the change in stream wise length scale during a pairing. The linear mechanism for this process was found to be very slow, requiring roughly three pairings to complete a doubling of the span wise scale. Stronger three-dimensionality can produce more rapid scale changes but is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one span wise scale to a similar array at a larger span wise scale